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Abstract. Algorithmic Cooling (AC) of Spins is potentially the first near-future application 
of quantum computing devices. Straightforward quantum algorithms combined with novel 
entropy manipulations can result in a method to improve the identification of molecules. 
We introduce here several new exhaustive cooling algorithms, such as the Tribonacci and 
fc-bonacci algorithms. In particular, we present the "all-bonacci" algorithm, which appears 
to reach the maximal degree of cooling obtainable by the optimal AC approach. 

1 Introduction 

Molecules are built from atoms, and the nucleus inside each atom has a property called "spin" . The 
spin can be understood as the orientation of the nucleus, and when put in a magnetic field, certain 
spins are binary, either up (ZERO) or down (ONE). Several such bits (inside a single molecule) 
represent a binary string, or a register. A macroscopic number of such registers/molecules can 
be manipulated in parallel, as is done, for instance, in Magnetic Resonance Imaging (MRI). The 
purposes of magnetic resonance methods include the identification of molecules (e.g., proteins), 
material analysis, and imaging, for chemical or biomedical applications. From the perspective of 
quantum computation, the spectrometric device that typically monitors and manipulates these 
bits/spins can be considered a simple "quantum computing" device. 

Enhancing the sensitivity of such methods is a Holy Grail in the area of Nuclear Magnetic 
Resonance (NMR). A common approach to this problem, known as "effective cooling", has been 
to reduce the entropy of spins. A spin with lower entropy is considered "cooler" and provides a 
better signal when used for identifying molecules. To date, effective cooling methods have been 
plagued by various limitations and feasibility problems. 

"Algorithmic Cooling" [11213] is a novel and unconventional effective-cooling method that 
vastly reduces spin entropy. AC makes use of "data compression" algorithms (that are run on 
the spins themselves) in combination with "thermalization" . Due to Shannon's entropy bound 
(source-coding bound [4]), data compression alone is highly limited in its ability to reduce en- 
tropy: the total entropy of the spins in a molecule is preserved, and therefore cooling one spin is 
done at the expense of heating others. Entropy reduction is boosted dramatically by taking ad- 
vantage of the phenomenon of thermalization, the natural return of a spins entropy to its thermal 
equilibrium value where any information encoded on the spin is erased. Our entropy manipulation 
steps are designed such that the excess entropy is always placed on pre-selected spins, called "reset 
bits" , which return very quickly to thermal equilibrium. Alternating data compression steps with 
thermalization of the reset spins thus reduces the total entropy of the spins in the system far be- 
yond Shannon's bound. The AC of short molecules is experimentally feasible in conventional NMR 
labs; we, for example, recently cooled spins of a three-bit quantum computer beyond Shannon's 
entropy bound [5] . 



1.1 Spin- Temperature and NMR Sensitivity 

For two-state systems (e.g. binary spins) there is a simple connection between temperature, en- 
tropy, and probability. The difference in probability between the two states is called the polarization 
bias. Consider a single spin particle in a constant magnetic field. At equilibrium with a thermal 



heat-bath the probabilities of this spin to be up or down (i.e., parallel or anti-parallel to the 
magnetic field) are given by: p-f = 1+ 2 e ° , and pi = 1 ~ £ ° . We refer to a spin as a bit, so that 
||) = |0) and \l) = |1), where \x) represents the spin-state x. The polarization bias is given by 

£ a = P\ ~ Pi = tanh ( ^k^t ) > wnere B i s the magnetic field, 7 is the particle-dependent gyromag- 
netic constantQ Kb is Boltzmann's coefficient, and T is the thermal heat-bath temperature. Let 
e = 2KbT sucn that £0 = tanhe. For high temperatures or small biases, higher powers of e can be 
neglected, so we approximate £0 ~ £■ Typical values of £0 f° r nuclear spins (at room temperature 
and a magnetic field of ~ 10 Tesla) are — 10~ 6 . 

A major challenge in the application of NMR techniques is to enhance sensitivity by over- 
coming difficulties related to the Signal-to-Noise Ratio (SNR). Five fundamental approaches were 
traditionally suggested for improving the SNR of NMR. Three straightforward approaches - cool- 
ing the entire system, increasing the magnetic field, and using a larger sample - are all expensive 
and limited in applicability, for instance they are incompatible with live samples. Furthermore, 
such approaches are often impractical due to sample or hardware limitations. A fourth approach 
- repeated sampling - is very feasible and is often employed in NMR experiments. However, an 
improvement of the SNR by a factor of M requires M 2 repetitions (each followed by a significant 
delay to allow relaxation), making this approach time-consuming and overly costly. Furthermore, 
it is inadequate for samples which evolve over the averaged time-scale, for slow-relaxing spins, or 
for non-Gaussian noise. 

1.2 Effective cooling of spins 

The fifth fundamental approach to the SNR problem consists of cooling the spins without cooling 
the environment, an approach known as "effective cooling" of the spins 6 7 8 9J. The effectively- 
cooled spins can be used for spectroscopy until they relax to thermal equilibrium. The following 
calculations are done to leading order in £0 and are appropriate for eo « 1. A spin temperature at 
equilibrium is T cx £q 1 . The single-spin Shannon entropy is H = 1— (sq/ In 4) . A spin temperature 
out of thermal equilibrium is similarly defined (see for instance [10]). Therefore, increasing the 
polarization bias of a spin beyond its equilibrium value is equivalent to cooling the spin (without 
cooling the system) and to decreasing its entropy. 

Several more recent approaches are based on the creation of very high polarizations, for exam- 
ple dynamic nuclear polarization [11] , para- hydrogen in two-spin systems [12] , and hyperpolarized 
xenon [13j . In addition, there are other spin-cooling methods, based on general unitary transfor- 
mations [7] and on (closely related) data compression methods in closed systems [8]. 

One method for effective cooling of spins, reversible polarization compression (RPC), is based 
on entropy manipulation techniques. RPC can be used to cool some spins while heating oth- 
ers |7I8] , Contrary to conventional data compression!! RPC techniques focus on the low-entropy 
spins, namely those that get colder during the entropy manipulation process. RPC, also termed 
"molecular-scale heat engine" , consists of reversible, in-place, lossless, adiabatic entropy manipu- 
lations in a closed system (8]. Therefore, RPC is limited by the second law of thermodynamics, 
which states that entropy in a closed system cannot decrease, as is also stated by Shannon's 
source coding theorem [4]. Consider the total entropy of n uncorrelated spins with equal biases, 
H{n) w n(l — £q/1u4). This entropy could be compressed into m > n(l — £ 2 /ln4) high entropy 
spins, leaving n — m extremely cold spins with entropy near zero. Due to preservation of entropy, 
the number of extremely cold spins, n — m, cannot exceed n£§/ln4. With a typical £0 ~ 10 -5 , 
extremely long molecules (~ 10 10 atoms) are required in order to cool a single spin to a tempera- 
ture near zero. If we use smaller molecules, with n <C 10 10 , and compress the entropy onto n — 1 
fully-random spins, the entropy of the remaining spin satisfies [2] 

1 This constant, 7, is thus responsible for the difference in the equilibrium polarization bias of different 
spins [e.g., a hydrogen nucleus is about 4 times more polarized than a 1 C nucleus, but less polarized 
by three orders of magnitude than an electron spin]. 

2 Compression of data [4] such as bits in a computer file, can be performed by condensation of entropy 
to a minimal number of high entropy bits, which are then used as a compressed file. 



1 - 4„ai > r»(l - ell In 4) - (n - 1) = 1 - ne§/ In 4. 



(1) 



Thus, the polarization bias of the cooled spin is bounded by 

Efinal < eoV™ . (2) 

When all operations are unitary, a stricter bound than imposed by entropy conservation was 
derived by S0rensen [7]. In practice, due to experimental limitations, such as the efficiency of the 
algorithm, relaxation times, and off-resonance effects, the obtained cooling is significantly below 
the bound given in eqO 

Another effective cooling method is known as polarization transfer fPT) ^o\T4\ . This technique 
may be applied if at thermal equilibrium the spins to be used for spectroscopy (the observed 
spins) are less polarized than nearby auxiliary spins. In this case, PT from the auxiliary spins to 
the observed spins is equivalent to cooling the observed spins (while heating the auxiliary spins). 
PT from one spin to another is limited as a cooling technique, because the polarization bias 
increase of the observed spin is bounded by the bias of the highly polarized spin. PT is regularly 
used in NMR spectroscopy, among nuclear spins on the same molecule [6]. As a simple example, 
consider the 3-bit molecule trichloroethylene (TCE) shown in Fig. [1] The hydrogen nucleus is 



CI CI 




Fig. 1. A 3-bit computer: a TCE molecule labeled with two 13 C. TCE has three spin nuclei: two 
carbons and a hydrogen, which are named CI, C2 and H. The chlorines have a very small signal 
and their coupling with the carbons is averaged out. Therefore, TCE acts as a three-bit computer. 
The proton can be used as a reset spin because relative to the carbons, its equilibrium bias is 
four times greater, and its thermalization time is much shorter. Based on the theoretical ideas 
presented in [2], the hydrogen of TCE was used to cool both carbons, decrease the total entropy 
of the molecule, and bypass Shannon's bound on cooling via RPC [5]. 

about four times more polarized than each of the carbon nuclei; PT from a hydrogen can be used 
to cool a single carbon by a factor of four. A different form of PT involves shifting entropy from 
nuclear spins to electron spins. This technique is still under development |9ll3j . but has significant 
potential in the future. 

Unfortunately, the manipulation of many spins, say n > 100, is a very difficult task, and the 
gain of y/n in polarization is not substantial enough to justify putting this technique into practice. 

In its most general form, RPC is applied to spins with different initial polarization biases, thus 
PT is a special case of RPC. We sometimes refer to both techniques and their combination as 
reversible algorithmic cooling. 

2 Algorithmic Cooling 

Boykin, Mor, Roychowdhury, Vatan, and Vrijen (hereinafter referred to as BMRVV), coined the 
term Algorithmic Cooling (AC) for their novel effective-cooling method pQ. AC expands previous 
effective-cooling techniques by exploiting entropy manipulations in open systems. It combines RPC 



with relaxation (namely, thermalization) of the hotter spins, in order to cool far beyond Shannon's 
entropy bound. 

AC employs slow- relaxing spins (which we call computation spins) and rapidly relaxing spins 
(reset spins), to cool the system by pumping entropy to the environment. Scheme 1 details the 
three basic operations of AC. The ratio R r elax—timesi between the spin-lattice relaxation times of 
the computation spins and the reset spins, must satisfy R r eiax-times 3* lj to permit the application 
of many cooling steps to the system. 

In all the algorithms presented below, we assume that the relaxation time of each computation 
spin is sufficiently large, so that the entire algorithm is completed before the computation spins 
lose their polarization. 

The practicable algorithmic cooling (PAC) suggested in [2] indicated a potential for near-future 
application to NMR spectroscopy [3]. In particular, it presented an algorithm (named PAC2) which 
uses any odd number of spins such that one of them is a reset spin and the other 2L spins are 
computation spins. PAC2 cools the spins such that the coldest one can (ideally) reach a bias of 
(3/2) L . This proves an exponential advantage of AC over the best possible reversible algorithmic 
cooling, as reversible cooling techniques (e.g., of refs [7] and [5J) are limited to a bias improvement 
factor of y/n. As PAC2 can be applied to small L (and small n) , it is potentially suitable for near 
future applications. 

Scheme 1: AC is based on the combination of three distinct operations: 

1. RPC. Reversible Polarization Compression steps redistribute the entropy in the system so 
that some computation spins are cooled while other computation spins become hotter than the 
environment. 

2. SWAP. Controlled interactions allow the hotter computation spins to adiabatically lose their 
entropy to a set of reset spins, via PT from the reset spins onto these specific computation 
spins. 

3. WAIT. The reset spins rapidly thermalize, conveying their entropy to the environment, while 
the computation spins remain colder, so that the entire system is cooled. 

2.1 Block- Wise Algorithmic Cooling 

The original Algorithmic Cooling (BMRVV AC) [T] was designed to address the scaling problem 
of NMR quantum computing. Thus, a significant number of spins are cooled to a desired level and 
arranged in a consecutive block, such that the entire block may be viewed as a register of cooled 
spins. The calculation of the cooling degree attained by the algorithm was based on the law of 
large numbers, yielding an exponential reduction in spin temperature. Relatively long molecules 
were required due to this statistical nature, in order to ensure the cooling of 20, 50, or more spins 
(the algorithm was not analyzed for a small number of spins). All computation spins were assumed 
to be arranged in a linear chain where each computation spin (e.g., 13 C) is attached to a reset 
spin (e.g., 1 H). Reset and computation spins are assumed to have the same bias £q. BMRVV AC 
consists of applying a recursive algorithm repeating (as many times as necessary) the sequence of 
RPC, SWAP (with reset spins), and WAIT, as detailed in Scheme 1. The relation between gates 
and NMR pulse sequences is discussed in ref [15] , 

This cooling algorithm employed a simple form of RPC termed Basic Compression Subroutine 
(BCS) [T]. The computation spins are ordered in pairs, and the following operations are applied 
to each pair of computation spins X, Y. Xj, Yj denote the state of the respective spin after stage 
j; X ,Yq indicate the initial state. 

1. Controlled-NOT (CNOT), with spin Y as the control and spin X as target: Y\ = Y ,Xi = 
Ao © Yq, where © denotes exclusive OR (namely, addition modulo 2). This means that the 
target spin is flipped (NOT gate ||) «-* \[)) if the control spin is |1) (i.e., \[)). Note that 
Xi = |0) Ao = Yq, in which case the bias of spin Y is doubled, namely the probability that 
Yi = |0) is very close to (1 + 2e )/2. 

2. NOT gate on spin X x . So X 2 = NOT(Xx). 



3. Controlled-SWAP (CSWAP), with spin X as a control: if X 2 = |1) (Y was cooled), transfer 
the improved state of Y (doubled bias) to a chosen location by alternate SWAP and CSWAP 
gates. 



Let us show that indeed in step 1 the bias of spin Y is doubled whenever X\ = |0). Consider the 
truth table for the CNOT operation below: 

input : X Y output : X{Y\ 
-> 
1 -> 11 

10 -> 10 

11 -> 1 

The probability that both spins are initially |0) is p oa = P (X = |0) , Y = |0» = (1 + £o) 2 / 4 - 
In general, p M = P (X = \k) ,Y = \l)) such that p 01 = p w = (1 + e )(l - £o)/4, and p n = 
(1 - £o ) 2 /4. After CNOT, the conditional probability P (ft = |0) |Xj = |0)) is = _J2aa_ 

where qki = P (Xi = \k) , ft — |Z)) and is derived from pj~i according to the truth table, so that 
<?oo = Poo j and qoi = pn. In terms of eo this probability is 

(1 + eo) 2 l + 2e 



(l + £o ) 2 + (l- £o ) 2 



indicating that the bias of Y was indeed doubled in this case. 

For further details regarding the basic compression subroutine and its application to many 
pairs of spins (in a recursive manner) to reach a bias of 2 J eo, we refer the reader to ref PQ. 



2.2 Practicable Algorithmic Cooling (PAC) 

An efficient and experimentally feasible AC technique was later presented, termed "practicable 
algorithmic cooling (PAC)" \2 . Unlike the first algorithm, the analysis of PAC does not rely on 
the law of large numbers, and no statistical requirement is invoked. PAC is thus a simple algorithm 
that may be conveniently analyzed, and which is already applicable for molecules containing very 
few spins. Therefore, PAC has already led to experimental implementations. PAC algorithms 
use PT steps, reset steps, and 3-bit- compression (3B-Comp). As already mentioned, one of the 
algorithms presented in [2], PAC2, cools the spins such that the coldest one can (ideally) reach a 
bias of (3/2) L , while the number of spins is only 2L + 1. PAC is simple, as all compressions are 
applied to three spins, often with identical biases. The algorithms we present in the next section 
are more efficient and lead to a better degree of cooling, but are also more complicated. We believe 
that PAC is the best candidate for near future applications of AC, such as improving the SNR of 
biomedical applications. 

PAC algorithms [2 use a basic 3-spin RPC step termed 3-bit-compression (3B-Comp): 
Scheme 2: 3-BitCompression (SBComp) 

1. CNOT, with spin B as a control and spin A as a target. Spin A is flipped if B = |1). 

2. NOT on spin A. 

3. CSWAP with spin A as a control and spins B and C as targets. B and C are swapped if 
A=\l). 

Assume that the initial bias of the spins is sq. The result of scheme 2 is that spin C is cooled: 
if A = |1) after the first step (and A = |0) after the second step), C is left unchanged (with its 
original bias eo); if however, A = |0) after the first step (hence A = |1) after the second step), spin 
B is cooled by a factor of about 2 (see previous subsection), and following the CSWAP the new 
bias is placed on C. Therefore, on average, C is cooled by a factor of 3/2. We do not care about 
the biases of the other two spins, as they subsequently undergo a reset operation. 



In many realistic cases the polarization bias of the reset spins at thermal equilibrium, £q, is 
higher than the biases of the computation spins. Thus, an initial PT from reset spins to compu- 
tation spins (e.g., from hydrogen to carbon or nitrogen), cools the computation spins to the th 
purification level, eq. 

As an alternative to Scheme 2, the 3B-Comp operation depicted in Scheme 3 is similar to the 
CNOT-CSWAP combination (Scheme 2) and cools spin C to the same degree. This gate is known 
as the MAJORITY gate since the resulting value of bit C indicates whether the majority of the 
bits had values of |0) or |1) prior to the operation of the gate. 
Scheme 3: Single operation implementing 3B-Comp 
Exchange the states 1 100) <-> 1 01 1) . 
Leave the rest of the states unchanged. 

If 3B-Comp is applied to three spins {C,B,A} which have identical biases, £c = £b = £a = £o, 
then spin C will acquire a new bias e' c . This new bias is obtained from the initial probability that 
spin C is |0) by adding the effect of exchanging 1 100) <-* 1 1 1) : 

1 + £ 'c 1 + ec . r „s 
— 2 — = — ^ — +2>|ioo) -P|ou> (3) 

1 + e 1 - e 1 + £q 1 + £o 1 + e 1 - £o 1 ~ 1 + 



The resulting bias is 

and in the case where £q C 1 



3eo — 
~ 2 

, 3£o 



(4) 



(5) 



We have reviewed two schemes for cooling spin C: one using CNOT and CSWAP gates (see 
Scheme 2), and the other using the MAJORITY gate (see Scheme 3). Spin C reaches the same 
final bias following both schemes. The other spins may obtain different biases, but this is irrelevant 
for our purpose, as they undergo a reset operation in the next stages of any cooling algorithm. 

The simplest practicable cooling algorithm, termed Practicable Algorithmic Cooling 1 (PAC1) [2], 
employs dedicated reset spins, which are not involved in compression steps. PAC1 on three com- 
putation spins is shown in Fig. [2j where each computation spin X has a neighboring reset spin, 
rx, with which it can be swapped. The following examples illustrate cooling by PAC1. In order to 
cool a single spin (say, spin C) to the first purification level, start with three computation spins, 
CBA, and perform the sequence presented in Example l.@ 
Example 1: Cooling spin C to the 1 st purification level by PAC1 

1. PT((rc — > C); (t~b — * B); (ta ^ A) ), to initiate all spins. 

2. 3B-Comp(C] B] A), increase the polarization ofC. 

A similar algorithm can also cool the entire molecule beyond Shannon's entropy bound. See 
the sequence presented in Example 2. 

Example 2: Cooling spin C and bypassing Shannon's entropy bound by PAC1 

1. PT({rc — > C); (tb — * B); (ta —> A)), to initiate all spins. 

2. 3B-Comp(C; B; A), increase the polarization ofC. 

3. WAIT 

Jf.. PT({tb — * B); (rA — * A) ), to reset spins A,B. 

In order to cool one spin (say, spin E) to the second purification level (polarization bias £2), 
start with five computation spins (EDCBA) and perform the sequence presented in Example 3. 

For small biases, the polarization of spin E following the sequence in Example 3 is £2 ~ 
(3/2) El w (3/2) 2 £o . 

Example 3: Cooling spin E in EDCBA to the 2nd purification level by PAC1 



3 This sequence may be implemented by the application of an appropriate sequence of radiofrequency 
pulses. 



1. PT(E; D;C), to initiate spins EDC. 

2. 3B-Comp(E; D;C), increase the polarization of E to £\. 

3. WAIT 

4- PT(D;C; B), to initiate spins DCB. 

5. 3B-Comp(D;C\ B), increase the polarization of D to E\. 

6. WAIT 

7. PT(C;B;A), to initiate spins CBA. 

8. 3B-Comp(C] B] A), increase the polarization of C to £\. 

9. 3B-Comp(E; D;C), increase the polarization of E toe?.. 

For molecules with more spins, a higher cooling level can be obtained. 

This simple practicable cooling algorithm (PAC1) is easily generalized to cool one spin to any 
purification level L. [5] The resultant bias will be very close to (3/2) L , as long as this bias is 
much smaller than 1. For a final bias that approaches 1 (close to a pure state), as required for 
conventional (non-ensemble) quantum computing, a more precise calculation is required. 

Consider an array of n computation spins, c„c n _i . . . C2C1, where each computation spin, Cj, is 
attached to a reset spin, n (see Fig. [2] for the case of n = 3). To cool Ck, the spin at index k, to a 
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Fig. 2. An abstract example of a molecule with three computation spins A, B and C, attached 
to reset spins tai r B, and r<7, respectively. All spins have the same equilibrium polarization bias, 
eo (a). The temperature after each step is illustrated by color: gray - thermal equilibrium; white - 
colder than initial temperature; and black - hotter than initial temperature. PAC uses 3-bit com- 
pression (3B-Comp), Polarization Transfer (PT) steps, and RESET steps: 1. 3B-Comp(C; B; A); 
the outcome of this step is shown in (b). 2. PT(rg — > B), PT(rc — > C); the outcome is shown 
in (c). 3. RESET (rs, rc)', the outcome is shown in (d). The 3-bit-compression applied in the first 
step operates on the three computation spins, increasing the bias of spin A by a factor of 3/2, 
while heating the other two spins. The 3B-Comp step cools spin A, and the following PT and 
RESET steps restore the initial biases of the other spins, thus the entire system is cooled. 



purification level 36 1...L the procedure Mj(k) was recursively defined as follows 0: M (k) is 
defined as a single PT step from reset spin rk to computation spin Ck to yield a polarization bias 
of eo (the th purification level). The procedure Mi(fc) applies Mq to the three spins followed by 
3B-Comp on these spins, so that spin Ck is cooled to the first purification level. Similarly, Mb (As) 
applies Mi three times to cool Ck', Ck—i', Ck—2 to the first purification level, followed by 3B-Comp 
on these spins, so that spin Ck is cooled to the second purification level. We use the notation 
B/j-u^j(k) to represent the application of 3B-Comp to spins to purify spin Ck from Ej-i to £j. 
Then, the full algorithm has a simple recursive form, described in Algorithm 1. 
Algorithm 1: Practicable algorithmic cooling 1 (PAC1): 
Forje{l,...,L} 

Mj(k) = 6 {0 _ 1) ^ } (fc)M J _ 1 (fc - 2) Mj-x{k - 1) Mj-i(k) , (6) 
applied from right to left (Mj-\{k) is applied first). 

For instance, Mi (3) = i3{ -»i}(3)Mo(l) Mo(2) Mq(3), is 3B-Comp applied after reset as described 



in Example 10 Clearly, M\(k) can be applied to any k > 3, M^k) to k > 5, and Mj(k) to 
k > 2j + 1. Thus, to cool a single spin to a purification level of L, 2L + 1 computation spins and 
an equal number of reset spins are required. A single reset spin could be used for initializing all 
relevant computation spins, at the expense of additional time steps. 

Reset spins may also be used for compression, thus replacing the 3B-Comp and PT steps 
above by a generalized RPC. The corresponding algorithm, termed PAC2, has an improved space 
complexity relative to PAC1. We explicitly show how this is achieved. Let £o be the polarization 
bias of the reset spin. In order to cool a single spin to Si, start with two computation spins, CB, 
and one reset spin, A, and perform the sequence shown in Example 4 to cool spin C . 
Example 4: Cooling spin C to the 1 st purification level by PAC2 

1. PT(A^ B). 

2. PT(B —* C) to initiate spin C. 

3. RESET(A) (by waiting). 

4- PT(A —tB)to initiate spin B. 

5. RESET(A). If the thermalization time of the computation spins is sufficiently large, there are 
now three spins with polarization bias Sq . 

6. 3B-Comp to increase the polarization of spin C to £\. 

In order to cool one spin (say, spin E) to the second purification level (polarization bias £2), 
start with 5 computation spins (EDCBA) and follow Example 5. 
Example 5: Cooling spin E in EDCBA to the 2 nd purification level by PAC2 

1. PT sequentially to initiate spins EDC (RESET(A) after each PT). 

2. 3B-Comp on spins EDC to increase the polarization of spin E to £\. 

3. PT sequentially to initiate spins DCB (RESET(A) after each PT). 
4- 3B-Comp on spins DCB to increase the polarization of spin D to E\. 

5. PT sequentially to initiate spins CB (RESET(A) after each PT). 

6. 3B-Comp on spins CBA to increase the polarization of spin C to £\. 

7. 3B-Comp on spins EDC to increase the polarization of spin E to £2. 

By repeated application of PAC2 in a recursive manner (as for PAC1), spin systems can be 
cooled to very low temperatures. PAC1 uses dedicated reset spins, while PAC2 also employs reset 
spins for compression. The simplest algorithmic cooling can thus be obtained with as few as 3 
spins, comprising 2 computation spins and one reset spin. 

The algorithms presented so far applied compression steps (3B-Comp) to three identical biases 
(£0); recall that this cools one spin to a new bias, e' c f=a (3/2)eo- Now consider applying compression 
to three spins with different biases (ec, £b, £a); spin C will acquire a new bias e' c , which is a 
function of the three initial biases |16ll7j . This new bias is obtained from the initial probability 
that spin C is |0) by adding the effect of the exchange 1 100) |011): 

l + e'c 1 + £ C . 

— = — 2 — hPl0 ° ~ Pou 

l + e c , 1 - £c 1 + £b 1 + £a 1 + £ C 1 - £b 1 - £a 



1 _i_ ec+EB+eA—eceseA 
1 ' 2 



The resulting bias is 

fc c - 

and in the case where sc, £s, £a *C 1 



2 

_ £c + £b + £a - s C £bSa 



(7) 

(8) 



e c + £b + £a ^ 



4 The procedure of cooling one spin to the second level (starting with five spins) is written as M?(5) = 
B { i^3}(5)Mi(3)ATi(4)Mi(5). 



3 Exhaustive Cooling Algorithms 



The following examples and derivations are to leading order in the biases. This is justified as long 
as all biases are much smaller than 1, including the final biases. For example, with Eq ~ 10~ 5 (sq is 
in the order of magnitude of 10 -5 ) and n < 13, the calculations are fine (see section|3]for details). 



3.1 Exhaustive Cooling on Three Spins 

Example 6: Fernandez 116]: T(C, B, A,m) 
Repeat the following m times 

1. 3B-Comp(C;B;A), to cool C . 

2. RESET(B-A) 

Consider an application of the primitive algorithm in Example 6 to three spins, C, -B, A, where 
A and B are reset spins with initial biases of and (the index over the bias of C denotes the 
iteration). After each WAIT step, spins A and B are reset back to their equilibrium bias, eo- As B 
and A play roles of both reset and computation spins, the operation RESET in Example 6 simply 
means WAIT. From eqEl after the first iteration 



_(i) _ 4 0) +£B+e A _ + 2e 



£o- 



After the second iteration 



After the m th iteration 



.(2) _ £ ( c +e B +e A _ £ + 2e = 3_£o 
~~ c 2 2 ~ 2 



(m— 1) n 

_( TO ) _ e c + 2e _ (o) 



2e J2 2 ~ J = + (1 - 2~ m ) 2e . (10) 



The asymptotic bias (m —> oo) may be extracted from 



(m) £ c ' ZS(] 



(11) 



with the unique solution and bias configuration 



4 m) = 2 £o ^{2£o,£o,£o}- (12) 

In order to achieve good asymptotics, one would like to reach ec = (2 — S)so, where S is arbitrarily 
small. In this case the number of iterations required is given by 2 1_m = 5 ==>• m = 1 + |~log 2 (1/ S)~\ . 
For example, if S = 10~ 5 , 18 repetitions are sufficient. Up to an accuracy of 8, the biases after the 
last reset are as in eg IT2"1 



3.2 The Fibonacci Algorithm 

An algorithm based on 3B-Comp was recently devised [18] , which produces a bias configuration 
that asymptotically approaches the Fibonacci series. In particular, when applied to n spins, the 
coldest spin attains the bias £fi na i ~ SoF n , where F n is the n th element of the series and SoF n <C 1. 
Note that of eq[T2]is £0^3- Also note that for 12 spins and £0 ~ 10 -5 , £0^12 *C 1, so the 
approximation is useful for non-trivial spin systems. Compare the bias enhancement factor in this 
case, F12 = 144, to PAC2 with 13 spins - (3/2) 6 » 11. 



Example 7 expands Example 6 to four spins with initial biases Ed — £ c = and Eb = £a = £o 
(A, B are reset spins). We refer to the parameter m from Example 6 as m.3, the number of 
repetitions applied to three spins. We refer to T from Example 6 as T2 for consistency with the 
rest of the section. Consider the following example: 
Example 7: J-2,{D, C, B, A, mi, ITI3) 
Repeat the following m\ times: 

1. 3B-Comp(D]C] B), places the new bias on spin D. 

2. F 2 (C,B,A,m 3 ). 

After each iteration, i, the bias configuration is {e^, £™ 3 , £0, £ o}- Running the two steps in 
Example 7 exhaustively (7774,7773 S> 1) yields: 

£ c ^ ~ 2 £ °' 

4 mi) » £B + £c - + £d (13) 

, A m i) _ - , -(Wis) _n. _ _ 771 
— Eb + £q — o£o — 

This unique solution was obtained by following the logic of eqs [11] and 1121 

We generalize this algorithm to n spins A n , . . . , A\ in Algorithm 2. 
Algorithm 2: Fibonacci 3-2{A n , . . . , A%, m n , , 7773) 
Repeat the following m n times: 

1. 3B-Comp(A n ;A n -i;A n ^ 2 )- 

2. J : 2{A n ^x,...,Ax,m n -i,---,m 5 ). 

[with ^-"2(^3, A2, A±, TO3) defined by Example 6.] 

Note that different values of m n _i, . . . ,7713 may be applied at each repetition of step 2 in Al- 
gorithm 2. This is a recursive algorithm; it calls itself with one less spin. Running Algorithm 2 
exhaustively (m„, m„_i, . . . , 777,3 ^> 1) results, similarly to eqQ2] in 



_( m i) 



2 

=*e£" 3 «e v . ( 14 ) 

This formula yields the Fibonacci series {. . . , 8, 5, 3, 2, 1, 1}, therefore ea { — » £o-Fi- We next devise 
generalized algorithms which achieve better cooling. An analysis of the time requirements of the 
Fibonacci cooling algorithm is provided in [18] . 



3.3 The Tribonacci Algorithm 

Consider 4-bit-compression (4B-Comp) which consists of an exchange between the states 1 1000) 
and 1 01 1 1) (the other states remain invariant similar to 3B-Comp). Application of 4B-Comp to 
four spins D,C, B, A with corresponding biases Eo^ sc, £b, £a 1 results in a spin with the 
probability of the state |0) given by 

— 2 — = — 2 — +P|iooo> -P|om> « 2 ' ' ^ 

following the logic of eqs [7] and [8] and finally 

e' D « (ea + eb + £c + 3e D )/4. (16) 

Example 8 applies an algorithm based on 4B-Comp to 4 spins, with initial biases Ed — £c = 
(A, £? are reset spins). In every iteration of Example 8, running step 2 exhaustively yields the 
following biases: Ec = 2eo,£_b = £a = £o- The compression step (step 1) is then applied onto the 
configuration 

Example 8: T^(D, C, B, A, 7714, 7773) 
Repeat the following 7774 times: 



1. 4B-Comp(D;C;B;A). 

2. F 2 (C,B,A,m 3 ). 

From eqUH e^ +1) = (4e + 3£^ } )/ 4 - For sufficiently large rrii and 7713 the algorithm produces 
final polarizations of 

(IH4) 1 ^ (1TI4) (>Tl4) A /T7\ 

£ d « £ o + ^eb £ !b » 4e o- (17) 

For more than 4 spins, a similar 4B-Comp based algorithm may be defined. Example 9 applies an 
algorithm based on 4B-Comp to 5 spins. 
Example 9: T 3 (E, D, C, B, A, m 5 , to 4 , m 3 ) 
Repeat the following 777,5 times: 

1. 4B-Comp(E;D;C\B). 

2. F 3 (D,C,B,A,m 4 ,m 3 ). 

Step 2 is run exhaustively in each iteration; the biases of DCBA after this step are £jj = 
4eo, £c = 2eo, B = A = eq. The 4B-Comp step is then applied to the biases of spins (E, D, C, B). 
Similarly to eq[T|)l £^ = (7£o + 3e^ )/4. Hence, for sufficiently large 7715 and 7774 the final bias 
of E is 

4 m5) « (7 £ o + 3e^)/4 => 4 m5) « 7 £o . (18) 
Algorithm 3 applies to an arbitrary number of spins n > 4. This is a recursive algorithm, which 
calls itself with one less spin. 

Algorithm 3: Tribonacci: J- 3 (A n , . . . , A\, m n , . . . , 777,3) 
Repeat the following m n times: 

1. 4B-Comp(A n ; A n ^i; A n - 2 ; A n - 3 ). 

2. ^(Ai-i, ■ ■ .,Ai,m n -i, ■ ■ ■ ,m 3 ). 

[With J- 3 {Ai, A 3l A2, A\, 777,4, m 3 ) given by Example 8. J 
The compression step, 4B-Comp, is applied to 



and results in 



(*) 

£ A n > £A„_d £ A„_ 2 > £ A„ 



= + £ ^- 2 + £ a„_ 3 + 3e£)/4. (19) 

For sufficiently large rrij , j = 3, . . . , n 



.W 



-3 ~r JC A„ 

£A„_! + £A„_ 2 + £A„_ 3 - (2°) 

For £^ 3 = 2e ,£A 2 = £ Ai = £0) the resulting bias will be £^™"' ) » £o^n; where T n is the 77 th 
Tribonacci number |f| As for the Fibonacci algorithm, we assume 6oT n <C 1. The resulting series is 
{...,24,13,7,4,2,1,1}. 



3.4 The fc-bonacci Algorithm 

A direct generalization of the Fibonacci and Tribonacci algorithms above is achieved by the ap- 
plication of (k + l)-bit-compression, or (k + l)B-Comp. This compression on k + 1 spins involves 
the exchange of 1 100 • • • 000) and |011 • ■ ■ 111), leaving the other states unchanged. When (k+ 1)B- 
Comp is applied to k + 1 spins with biases {eA k+1 , £A fc , • ■ ■ , £a 25 e Ai} 5 where A\ and A 2 are reset 
spins, the probability that the leftmost spin is |0) becomes (similarly to eg [T5)) 



(2 fc ' 1 -i)^, + _ 1 +E J t =1 ^ J 
+ Pioo---ooo — Poii-'-ni ~ 3 ■ (21) 



1 + _ l + £fe+i , _ 1 + 



5 The Tribonacci series (also known as the Fibonacci 3-step series) is generated by the recursive formula 
Oj = Qi-i + Qi-2 + cti— 3, where 03 = 2 and a% = tti = 1. 



Therefore, the bias of the leftmost spin becomes 

, „ (2 fc - 1 -l)^ fc+1 +E 3 t l£ A 3 

£ k+l jjfcZi ' \ ZZ ) 

Example 10 applies an algorithm based on [k + l)B-Comp to k + 1 spins. T k on k + 1 spins calls 
(recursively) — 1 on k spins (recall that F 3 on four spins called F 2 on three spins). 
Example 10: T k (A k+1 ,...,A 1: m k+1 , . , . , m 3 ) 
repeat the following rrik+i times: 

1. (k + l)B-Comp(A k+1 ,A k , A 2 ,A 1 ) 

2. ^ 7 / s _i(Afc,...,Ai,mfe,...,m 3 ) 

If step [2] is run exhaustively at each iteration, the resulting biases are 

£A k = 2 fe 2 £oj £ A k _ 1 — 2 fc 3 £o, ■ ■ ■ , £a 3 = 2eo, £a 2 — e A t = £o- 

The (fc + l)B-Comp step is then applied to the biases , 2 fe ~ 2 £o, 2 fe_3 £o, ■ ■ ■ , 2eo> £0j £o- From 

eq[22j 

.(ifD^C 2 *- 1 - 1 )^^ 2 *- 1 ^ 



Hence, for sufficiently large m^, J = 3, . . . , k + 1, the final bias of A k+ \ is 

(2 fc - i ~i)^:ir ) +2 fe - i £ o 



« => e^ tfc+u « 2*- 1 e . (23) 



For more than k + 1 spins a similar (fc + l)B-Comp based algorithm may be defined. Example 11 
applies such an algorithm to k + 2 spins. 
Example 11: F fc (A fc+ 2, . . . ,Ai,m k+2 , ■ ■ . ,m 3 ) 
Repeat the following steps m k + 2 times: 

1. {k + l)£-CWp(A fc+ 2, A fe+ i, . . . , A 3 , 

2. JFfc (-Afc+i , . . . , Ai,m k+ i, . . . ,7713) [defined in Example 11.] 

When step [5] is run exhaustively at each iteration, the resulting biases are 

EA fc+ i = 2 fe 1 Eq, £A k = 2 fe 2 £o, ■ ■ ■ ) £A 3 = %£0i £ A 2 = £ Ai = £q- 



The (k + l)B-Comp is applied to the biases. From Eq 

.(m) _ - ^% +2 + (2 fc - 1) 



£ A k+2 2 fe -i 
Hence, for sufficiently large rrij,j = 3, . . . , k + 2, the final bias of A k+2 is 

, v (2 k ~ 1 - l)e { T k+2) + (2 fe - l)e , , 
e { Z k + f « Ak 2 + k \ =► e { Zlf « (2 fc - 1) E0 . (24) 

Algorithm 4 generalizes Examples 10 and 11. 
Algorithm 4: k-bonacci: T k {A n , . . . , Ai,m n , . . . , m 3 ) 
Repeat the following m n times: 

1. {k + l)B-Comp(A n ,A n _ l , . . .,A n _ k ). 

2. ffc(i n -i,..., J 4i,m„_i,...,m 3 ). 



[with TkiAk+i, . . . , j4i,mfc_|_i, . . . ,7713) defined in Example 10.] 

The algorithm is recursive; it calls itself with one less spin. The compression step, (k + l)B-Comp, 
is applied to 

0) 

£ A n ' £ A n -n ■ ■ ■ £ A n - k ■ 

From Eq. [2H the compression results in 

(25) 



For sufficiently large m n , m n _i, etc. 



(2 fe - 1 -i) £ £" ) +E?= 1 ^„_ j 



2 fc-i 

fe 



This set of biases corresponds to the fc-step Fibonacci sequence which is generated by a recursive 
formula. 

ai,a 2 = l,a^ = < > . (27) 

lEi=i^-i ^>^ + l J 

Notice that for 3 < £ < k + 1, 

e-i 

a f = ^ ai = l + l + 2 + 4 + -- - + 2 £ - 4 + 2 e ~ 3 = 2 e ~ 2 . (28) 

i=l 

The algorithm uses £-bit-compression (/?-B-Comp) gates, where 3 < I < k + 1. 



3.5 The All-bonacci Algorithm 

In Example 11 (Tk applied to k+1 spins), the resulting biases of the computation spins, 2 k ~ 1 Eo, 2 k ~ 2 Eo, . . . , 2eq, Eq, 
were proportional to the exponential series {2fc_i,2fc_2, ... ,4, 2, 1}. For example, Ti on 3 spins 
results in {2£ ,£o,£o} (see Example 6), and T-$ on 4 spins results in {4e , 2e , So, £0} (see Ex- 
ample 8). This coincides with the fc-step Fibonacci sequence, where ag = 2 e ~ 2 , for a\ = 1 and 
I = 2, 3, k + 1 (see eg [28]) . This property leads to cooling of n + 1 spins by a special case of 
fc-bonacci (Algorithm 5), where fc-bonacci is applied to k = n — 1 spins. 
Algorithm 5: All-bonacci: TM\{A n , ■ • • , A±, m n . . . , 7713) 
Apply T n -i{A n , . . . , Ai,m n , ■ ■ ■ ,m 3 ). 

The final biases after all-bonacci are ea^ — » £o2 1-2 for i > 1. The resulting series is {. . . , 16, 8, 4, 2, 1, 1}. 

The all-bonacci algorithm potentially constitutes an optimal AC scheme, as explained in sec- 
tion n 



3.6 Density Matrices of the Cooled Spin Systems 

For a system of three spins in a completely mixed state (CMS), the density matrix of each spin is 
For a spin system in a completely mixed state, the density matrix of each spin is: 



^■J(g)X(g)X, where 



2 2 V 1 , 

The density matrix of the entire system, which is diagonal, is given by the tensor product pcm = 

Diag(pcM) = 2- 3 (l, 1, 1, 1, 1, 1, 1, 1). (30) 



For three spins in a thermal state, the density matrix of each spin is 

) • (31) 

and the density matrix of the entire system is given by the tensor product p? — ® p*p ® Pr ' 
This matrix is also diagonal. We write the diagonal elements to leading order in £o : 

Diag(pT) = 2 _3 (1 + 3e , 1 + £o, 1 + £o, 1 + eo, 1 - Eo, 1 - so, 1 - £q, 1 - 3s ). (32) 

Consider now the density matrix after shifting and scaling, p' — 2 n (p — 2~ n T)/eo. For any 
Diag{p) = (p 1 ,p 2 ,p 3 , ■ ■ ■ ,p n ) the resulting diagonal is Diag(p') = (pi,p' 2 ,Pa, . . ■ ,p' n ) with p'j = 
2 n (pj — 2~™)/eo- The diagonal of a shifted and scaled (S&S) matrix for a completely mixed state 
(eq|50D is 

Diag(p' CM ) = (0,0,0,0,0,0,0,0), (33) 

and for a thermal state (eq. [32]) 

Diag(p' T ) = (3, 1, 1, 1, -1, -1, -1, -3). (34) 

In the following discussion we assume that any element, p, satisfies p'eo -C 1. When applied to 
a diagonal matrix with elements of the form p = ^ (1 ztp'^o), this transformation yields the 
corresponding S&S elements, ±p'i. 

Consider now the application of the Fibonacci to three spins. The resultant bias configuration, 
{2eo,£o:£o} feq I12[) is associated with the density matrix 

pa4( 1+2 v 2 >( 1+ v>( 1+ v0- <-) 

with a corresponding diagonal, to leading order in eo, 

Diag (p$ 6 ) = 2- 3 (1 + 4e , 1 + 2e , 1 + 2e , 1,1,1- 2e , 1 - 2e , 1 - 4e ) . (36) 
The S&S form of this diagonal is 

Diag (p'^ b ) - (4, 2, 2, 0, 0, -2, -2, -4) . (37) 
Similarly, the S&S diagonal for Tribonacci on four spins is 

Diag (p' { T ] nb ) = (8, 6, 6, 4, 4, 2, 2, 0, 0, -2, -2, -4, -4, -6, -6, -8) . (38) 
and the S&S form of the diagonal for all-bonacci on n spins is 

Diag (pfi) = [2"-\ (2"" 1 - 2), (2"" 1 - 2), . . . , 2, 2, 0, 0, . . . , -2" - l] . (39) 
which are good approximations as long as 2™£o <C 1. 




Partner Pairing Algorithm Recently a cooling algorithm was devised that achieves a superior 
bias than previous AC algorithms. [18119] This algorithm, termed the Partner Pairing Algorithm 
(PPA), was shown to produce the highest possible bias for an arbitrary number of spins after any 
number of reset steps. Let us assume that the reset spin is the least significant bit (the rightmost 
spin in the tensor-product density matrix). The PPA on n spins is given in Algorithm 6. 
Algorithm 6: Partner Pairing Algorithm (PPA) 
Repeat the following, until cooling arbitrarily close to the limit. 



1. RESET - applied only to a single reset spin. 



2. SORT - A permutation that sorts the 2 n diagonal elements of the density matrix by decreasing 
value, such that po is the largest, andp2"-i is the smallest. 



Written in terms of e, the reset step has the effect of changing the traced density matrix of the 
reset spin to 



f 



1/1 + EO 



1 - e 



(40) 



for any previous state. From eq20]it is clear that Eq = tanhe as stated in the introduction. For a 
single spin in any diagonal mixed state: 



Po 



Pi 



RESET Po+Pl ( 1 + £o 



1 - E 



1 ( 1 + 60 



1 - e 



(41) 



For two spins in any diagonal state a reset of the least significant bit results in 



Po 



\ 



Pi 



P2 



RESET 



PO+Pl 



P3 J 

V 



^(l-eo) 



^(1 + e ) 



(42) 



*4**(l-so)/ 



Algorithm 6 may be highly inefficient in terms of logic gates. Each SORT could require an exponen- 
tial number of gates. Furthermore, even calculation of the required gates might be exponentially 
hard. 

We refer only to diagonal matrices and the diagonal elements of the matrices (applications of 
the gates considered here to a diagonal density matrix do not produce off-diagonal elements) . For 
a many-spin system, RESET of the reset spin, transforms the diagonal of any diagonal density 
matrix, p, as follows: 



Diag{p) = (po,pi,P2,P3, ...)-► 

Po+Pi ri . s Po+ft n , P2+R,, , P2+P3 n v 



(43) 



as the density matrix of each pair, pi and Pi+i (for even i) is transformed by the RESET step as 
described by eq|40]above. We use the definition of S&S probabilities, p' = 2 n (p — 2— n)/e - The 
resulting S&S diagonal is 



Diag(p') 



Po+Pi 
2 

Po+Pi 



+ (l+e )-2- 
2 „ |+2 nPo+Pi 



2" 
so 
2" 

£o 

Pi + P'l , 1 PO + Pi , 2 nP0±Pl PO + Pi _ 2 n P0+Pl 

2 '2 2 ' 2 2 '' 



(44) 



where the second element is shown in the final expression. We now use p — 2—n{sop' + 1), to 
obtain 



2 nP*±Pi = tl 2-"( £o p[, + l) + 2-"( e oK + l) 

2 2 

£o(Po+p[) + 2 p' +pf 1 
= ^ — 1 + £ o 7; • 



(45) 
(46) 



Ref [18] provides an analysis of the PPA. We continue, as in the previous subsection, to analyze 
the case of p'e *C 1, which is of practical interest. In this case 1 + eo Po 1 2 Pl ~ 1. Hence, the effect 
of a RESET step on the S&S diagonal is: 

Diag(p') = (p' ,p , 1 ,p , 2 ,p' 3 ,...)^ (47) 

Po + Pi , 1 Po + Pi _ , p'i + Pz , , P2 + Pa _ -, 
2 + ' 2 ' 2 + ' 2 

Consider now three spins, such that the one at the right is a reset spin. Following ref [T5], we 
initially apply the PPA to three spins which are initially at the completely mixed state. The first 
step of the PPA, RESET feci |47|) . is applied to the diagonal of the completely mixed state (cq [33)) . 
to yield 

Diag(p' CMS ) -> Diag(p' RESET ) = (1, -1, 1,-1, 1, -1, 1, -1). (48) 
This diagonal corresponds to the density matrix 

PRESET = ( 1 + £ ° X _ £ J , (49) 

namely to the bias configuration {0, 0, Eq}. The next PPA step, SORT, sorts the diagonal elements 
in decreasing order: 

Diag(p' SORT ) = (1, 1, 1, 1,-1, -1, -1, -1), (50) 
that arises from the density matrix 

PSORT = ^( 1 + £ \ B (51) 



2 3 V 1-eo. 

which corresponds to the biases {eo,0, 0}. The bias was thus transferred to the leftmost spin. In 
the course of repeated alternation between the two steps of the PPA, this diagonal will further 
evolve as detailed in Example 12. 

The rightmost column of Example 12 lists the resulting bias configurations. Notice that after 
the 5 th step the biases are identical. Also notice that in the 6 th and 10 th steps the states 1 100) and 
1 1 1) are exchanged (3B-Comp); after both these steps, the state of the system cannot be written 
as a tensor product^ Also notice that after step 13, the PPA applies a SORT which will also 
switch between these two states. The PPA, applied to the bias configuration {te 7£o,£o}, where 
1 < t < 2 is simply an exchange 1 100) <-* 1 01 1) . This is evident from the diagonal 

Diag(p') = (f + 2, t, t, t - 2, -t + 2, -t, -t, -t - 2). (52) 

Thus, the PPA on three spins is identical to the Fibonacci algorithm applied to three spins (see 
Example 6). This analogy may be taken further to some extent. When the PPA is applied onto four 
spins, the outcome, but not the steps, is identical to the result obtained by Tribonacci algorithm 
(see Example 8). These identical results may be generalized to the PPA and all-bonacci applied 
onto n > 3 spins. 

Example 12: Application of PPA to 3 spins 



This is due to classical correlations (not involving entanglement) between the spins; an individual bias 
may still be associated with each spin by tracing out the others, but such a bias cannot be interpreted 
as temperature. 







Diag(p') 


= (0,0,0,0,0,0,0,0) 


{0,0,0} 


step 


1 


RESET 


/ 1 11 11 11 1 \ 

(1,-1,1,-1, 1,-1,1, -1J 


{0,0,£ } 


step 


2 


SORT 


(1,1,1,1,-1,-1,-1,-1) 


{eo,0,0} 


step 


Q 



RESET 
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/o r\ o n n on o^ 
(2,0,2,0,0,-2,0,-2) 


f _ A _ 1 

{e , 0, £ } 


step 


4 


SORT 
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(2,2,0,0,0,0,-2,-2) 


{eo,£o,0} 


step 


5 


RESET 
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/Oil 11 1 i o\ 

(3,1,1,-1,1,-1,-1,-3) 


{£o,£o,£or 


step 
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SORT 
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(3,1,1,1,-1,-1,-1,-3) 


A. A. 


step 
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RESET 

■ — ■ > 


^O, 1, z, u, u, — z, — 1, — oj 


r 3e eo jr „l 


step 


8 


SORT 
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(3,2,1,0,0,-1,-2,-3) 


r 3e c eo 1 


step 


9 


RESET 


1(7,3,3,-1,1,-3,-3,-7) 


{^,£0,^0} 


step 


10 


SORT 
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1(7,3, 3, 1,-1,-3,-3,-7) 


N.A. 


step 


11 


RESET 
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±(7,3,4,0,0,-4,-3,-7) 




step 


12 


SORT 
> 


±(7,4,3,0,0,-3,-4,-7) 


{ 7 -f ,£o,^f} 


step 


13 


RESET 


±(15,7,7,-1,1,-7,-7,-15) 


{^,£0^0} 



4 Optimal Algorithmic Cooling 

4.1 Lower Limits of the PPA and All-bonacci 

Consider an application of the PPA to a three-spin system at a state with the S&S diagonal of 
eql37l This diagonal is both sorted and invariant to RESET, hence it is invariant to the PPA. Eql37l 
corresponds to the bias configuration {2eo,£o,£o} which is the limit of the Fibonacci algorithm 
presented above (eq [T2|) . This configuration is the "lower" limit of the PPA with three spins in 
the following sense: any "hotter configuration" is not invariant and continues to cool down during 
exhaustive PPA until reaching or bypassing this configuration. For four spins the diagonal of Eg 1381 
is invariant to the PPA for the same reasons. This diagonal corresponds to the bias configuration 
{4eo, 2eo, £o, £o}, which is the limit of the Tribonacci algorithm (or the all-bonacci) applied to four 
spins. Any "hotter configuration" is not invariant and cools further during exhaustive PPA, until 
it reaches (or bypasses) this configuration. 

Now, we follow the approximation of very small biases for n spins, where 2™£o -C 1. The 
diagonal in eq[39]is invariant to the PPA. This diagonal corresponds to the bias configuration 

{2™ 2 £o,2™ 3 £o,2" 4 £o, ■ ■ ■ , 2eo, £o, £o}, (53) 

which is the limit of the all-bonacci algorithm. As before, any "hotter configuration" is not invariant 
and cools further. It is thus proven that the PPA reaches at least the same biases as all-bonacci. 

We conclude that under the assumption 2 n £o 1, the n-spin system can be cooled to the bias 
configuration shown in eq 1531 When this assumption is not valid (e.g., for larger n with the same 
£o), pure qubits can be extracted; see theorem 2 in [19] (theorem 3 in ref [18]). 

Yet, the PPA potentially yields "colder configurations". We obtained numerical results for 
small spin systems which indicate that the limits of the PPA and all-bonacci are identical. Still, 
since the PPA may provide better cooling, it is important to put an "upper" limit on its cooling 
capacity. 

4.2 Upper Limits on Algorithmic Cooling 

Theoretical limits for cooling with algorithmic cooling devices have recently been established [18119] . 
For any number of reset steps, the PPA has been shown to be optimal in terms of entropy ex- 
traction (see ref [19] and more details in section 1 of ref 18 ). An upper bound on the degree of 
cooling attainable by the PPA is therefore also an upper bound of any AC algorithm. The following 
theorem regards a bound on AC which is the loose bound of ref [18] . 



Theorem 1 No algorithmic cooling method can increase the probability of any basis state to abov ^\ 
min{2~"e 2 e , 1}, where the initial configuration is the completely mixed stated This includes the 
idealization where an unbounded number of reset and logic steps can be applied without error or 
decoherence. 

The proof of Theorem Q] involves applying the PPA and showing that the probability of any state 
never exceeds 2~™e 2 £ . 

5 Algorithmic Cooling and NMR Quantum Computing 

We have been using the language of classical bits, however spins are quantum systems; thus, spin 
particles (two- level systems) should be regarded as quantum bits (qubits). A molecule with n spin 
nuclei can represent an n-qubit computing device. The quantum computing device in this case is 
actually an ensemble of many such molecules. In ensemble NMR quantum computing 20 21 22J 
each computer is represented by a single molecule, such as the TCE molecule of Fig. [T] and the 
qubits of the computer are represented by nuclear spins. The macroscopic number of identical 
molecules available in a bulk system is equivalent to many processing units which perform the 
same computation in parallel. The molecular ensemble is placed in a constant magnetic field, 
so that a small majority of the spins are aligned with the direction of the field. To perform a 
desired computation, the same sequence of external pulses is applied to all molecules/computers. 
Any quantum logic-gate can be implemented in NMR by a sequence of radio-frequency pulses 
and intermittent delay periods during which spins evolve under coupling [23] . Finally, the state 
of a particular qubit is measured by summing over all computers/molecules. The process of AC 
constitutes a simple quantum computing algorithm. However, unlike other quantum algorithms, 
the use of quantum logic gates does not produce any computational speed-up, but instead generates 
colder spins. This constitutes the first near- future application of quantum computing devices. AC 
may also have an important long-term application; it may enable quantum computing devices 
capable of running important quantum algorithms, such as the factorization of large numbers [24] . 
NMR quantum computers 20 21 22J are currently the most successful quantum computing devices 
(see for instance ref [25]), but are known to suffer from severe scalability problems [1126127] , AC 
can be used for building scalable NMR quantum computers of 20-50 quantum bits if electron spins 
are used for the PT and RESET steps. PT with electron spins |9|11] can enhance the polarization 
by three or four orders of magnitude. Unfortunately, severe technical difficulties have thus far 
impeded common practice of this technique. An example of such a difficulty is the need to master 
two very different electromagnetic frequencies within a single machine. However, in case such PT 
steps come into practice (using machinery that allows conventional NMR techniques as well), 
AC could be applied with much better parameters; First, £o could be increased to around 0.01- 
0.1. Second, the ratio Rrelax-times could reach 10 3 — 10 4 . With these figures, scalable quantum 
computers of 20-50 qubits may become feasible. 

6 Discussion 

Algorithmic Cooling (AC) harnesses the environment to enhance spin polarization much beyond 
the limits of reversible polarization compression (RPC). Both cooling methods may be regarded as 
a set of logic gates, such as NOT or SWAP, which are applied onto the spins. Polarization transfer 
(PT), for instance, a form of RPC, may be obtained by a SWAP gate. AC algorithms are composed 
of two types of steps: reversible AC steps (RPC) applied to two spins or more, and reset steps, in 
which entropy is shifted to the environment through reset spins. These reset spins thermalize much 
faster than computation spins (which are cooled), allowing a form of a molecular heat pump. A 
prerequisite for AC is thus the mutual existence (on the same molecule) of two types of spins with 

7 A tighter bound, po < 2~™e 2 e , was claimed by theorem 1 of ref [18] . 

8 It is assumed that the computation spins are initialized by the polarization of the reset spins. 



a substantial relaxation times ratio (of at least an order of magnitude). While this demand limits 
the applicability of AC, it is often met in practice (e.g., by 1 H vs 13 C in certain organic molecules) 
and may be induced by the addition of a paramagnetic reagent. The attractive possibility of using 
rapidly-thcrmalizing electron spins as reset bits is gradually becoming a relevant option for the 
far future. 

We have surveyed previous cooling algorithms and suggested novel algorithms for exhaustive 
AC: the Tribonacci, the /c-bonacci, and the all-bonacci algorithms. We conjectured the optimality 
of all-bonacci, as it appears to yield the same cooling level as the PPA of refs [18119) . Improving the 
SNR of NMR by AC potentially constitutes the first short-term application of quantum computing 
devices. AC is further accommodated in quantum computing schemes (NMR-based or others) 
relating to the more distant future [1 28 29 30J. 
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